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Abstract 

We compute the partial width of the FCNC top quark decay t — > ch in the 
framework of the Minimal Super symmetric Standard Model, where h = h°, H°, A is 
any of the neutral Higgs of the MSSM. We include the SUSY electroweak, Higgs, and 
SUSY-QCD contributions. Our results substantially improve previous estimations 
on the subject, and we find that there is a possibility that they can be measured at 
LHC. 



^alk presented at the International Workshop on quantum effects in the MSSM, Barcelona, 9-13 
September, 1997. To appear in the proceedings, World Scientific, ed. J. Sola. 
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1 Introduction 



We perform the computation of the Flavour Changing Neutral Current (FCNC) decay of 
the top quark into a charm quark and a neutral Higgs particle in the framework of the 
Minimal Supersymmetric Standard Model (MSSM), t — > ch where h is any of the neutral 
Higgs particles of the MSSM, namely the light or heavy CP-even ("scalar") Higgs bosons 
(h°, H°), or the CP-odd ( "pseudoscalar" ) Higgs (^4°)0- We compute the contributions 
from the SUSY electroweak, Higgs, and SUSY-QCD sectors, in a sparticle mass model 
motivated by model building and Renormalization Group Equations (RGE). However, 
we neither restrict ourselves to a spectrum of any SUSY-GUT model (such as SUGRA) 
-which would constrain the masses in a narrow range-, nor to a generic, phenomenological 
motivated, spectrum -which would have too many parameters to play with. 

There exist some computations of FCNC top quark decays, both in the SM and in 
the MSSM0, |], [5]. There has been some work concerning the decay channel into gauge 
bosons (t — > cV, V = 7, Z, g), see for example Ref.|§] for some works on the subject. The 
conclusion of these works is that the branching ratio of this decay is at most 10~ 5 , maybe 
a bit larger in the gluon channel. However, to our knowledge, there are not so many 
works on the FCNC decay of the top quark into Higgs in the MSSM||, and they are not 
so complete as in the case of the gauge bosons. For example in p it is concluded that 
the branching ratio for the decay channel t — > ch in the MSSM is at most of 10~ 9 , for the 
SUSY electroweak contributions, and 10~ 5 , for the SUSY-QCD contributions. However 
we think that the work of is not complete. They do not include effects of the Higgs 
particles in the loops, and they do not take into account the qi <1r h vertices, so they 
miss the potentially large contributions coming from the trilinear soft-SUSY-breaking 
terms At,b, and from the higgsino mass parameter /i. We find that a full treatment of 
the SUSY-QCD contributions may greatly enhance the FCNC width by some orders of 
magnitude. Therefore, a more general, and more rigorous, computation of the decay 
t — > ch is mandatory 0. 

In section ^ we make a summary of the technics of the computation, in sections [S] 
and|] we present our results for the SUSY electroweak and the SUSY-QCD contributions 
to the decay width t — > ch respectively, and finally we present the conclusions. 

2 Relevant lagrangians and form factors 

The computation of FCNC processes at one loop, unlike other calculations, does not 
involve renormalization of parameters, or wave functions, so one is left only with the 
computation of the different diagrams that contribute to the process. The generic one 
loop diagrams contributing to the decay under study are in fig. |I[ The vertex diagram 
V follows after a straightforward calculation. As for the diagrams St and S c we define a 
mixed self-energy, 

E tc (A;) = ft Y, L {k 2 ) P L + ft Y, R {k 2 ) P R + m t (s Ls (fc 2 ) P L + Y, Rs {k 2 ) P R ) (1) 

-where the m t factor multiplying the scalar part is arbitrary, put there only to maintain 
the same units between the different and compute the effects of Sj to the amplitude. 

2 We use a widely known notation[jl); see Refs.Q for further details and conventions. 
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We compute the different contributions to V, S t and S c and define an "effective" vertex 



-i T = -igu c {p) (F L P L + F R P R ) u t (k) . (2) 

We have taken into account all three generations of quarks and squarks, and have 
performed the usual checks of the computation, in particular that the form factors F^ and 
Fft are free of divergences before adding up the three quark generations, both analytically 
and numerically in the implementation of the code. 

After squaring the matrix element (||), and multiplying by the phase space factor, we 
can compute the decay width, and define the ratio 

B^c^J^L (3) 

which will be the main object under study. This ratio is not the total branching fraction 
of this decay mode, as there are many other channels that should be added up to the 
denominator of (f|) in the MSSM, such as the two and three body decays of the top quark 
into SUSY particles, and also the decay channel t — > H + 60, |7|]. For the mass spectrum 
used in the numerical analysis (see sections ||| and f|) the former are phase space closed, 
whereas the latter could have a large branching ratio. 

We would like to single out two pieces of the interaction lagrangian, namely the ones 
involving higgsino-sbottom-charm and Higgs-bottom-charm 

£~hbc = -gV cb c(R la \ c P L + R 2a hPR)X + ba 

£nbc = t=— V cb c(m c cot (3 P L + m b tan (3 P R ) bH + , (4) 
V 2M\y 

where V cb is the CKM matrix element, and A C; b are the charm and bottom Yukawa cou- 
plings^. Looking at these two pieces of the lagrangian one can get an estimation of the 
relative importance of the different form factors (|2]) (see section [3]). 



3 SUSY-EW contributions 

For the electroweak contributions to the decay channel t — > c h we work in the so called 
Super-CKM basis, that is, we take the simplification that the squark mass matrix diag- 
onalizes as the quark mass matrix, so that FCNC processes appear at one loop through 
the charged sector (charged Higgs and charginos) with the same mixing matrix elements 
as in the Standard Model (the CKM matrix). We introduce, as usual, the left-right mass 
matrix mixing elements between squarks0: m q (A q — /xjtan/?, cot /?}). 

We have taken into account the contributions from charginos (xt) an d down type 
squarks (d a , a = 1, 2, . . . , 6 = di, d 2 , . . . , b 2 , the mass eigenstates down squarks), and 
from charged Higgs and Goldstone bosons (H + ,G + ) and down type quarks (d,s,b). We 
have not included the diagrams with gauge bosons {W + ) as the largest contributions will 
come from the Yukawa couplings of the top and (at large tan 0) bottom quarks. However, 
the leading terms from longitudinal W + are included through the inclusion of Goldstone 
bosons. 

3 We refer again to j| for conventions and notation. 
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The input parameters chosen to illustrate the results in figs. are: 



tan (5 


= 35 


/' 


= -100 GeV 


M 


= 150 GeV 


M A o 


= 60 GeV 




= 150 GeV 


m~ bl = m~ q 


= 200 GeV 


A t = A q 


= 300 GeV 


A b 


= -300 GeV 



(5) 



where , m~ bi are the lightest t and b mass, and all the masses are above present experi- 
mental bounds ||-[[10], though the mass of the pseudoscalar Higgs (M40) is in the edge of 
present LEP bounds||, and it is possibly excluded by the last analyses on the charged 
Higgs mass]?], [L(| [11], However, this light Higgs mass is not essential in the results, as 
can bee seen in fig. ||] (d). We have chosen a SUSY mass spectrum around 200 GeV, which 
is not too light, so the results will not be artificially optimized. We have also checked all 
through the numerical analysis that other bounds on experimental parameters (such as 
Sp) are fulfilled. 

In fig. |2| we have plotted the different form factors of (fj) as a function of tan f3 for the 
channel with the lightest scalar Higgs (h°). We can see that the contributions from the 
Higgs sector and the contributions from the chargino sector are of the same order. It turns 
out that they can be either of the same sign, or of opposite sign. The chosen negative 
value for Ab is to make the two contributions of the same sign. It is also clear that in both 
cases Fr^> Fl. This can be easily understood by looking at the interaction lagrangians 
involving higgsino-sbottom-charm and Higgs-bottom-charm (^) where we can see that in 
both of them the contribution to the right-handed form factor will be enhanced by the 
Yukawa coupling of the bottom quark, compared with the charm Yukawa coupling that 
will contribute to the left-handed form factor. On the other hand we have checked that 
the inclusion of the first two generations of quarks and squarks only has an effect of a few 
percent on the total result. 

In fig. |3] we can see the evolution of the ratio (Q) with various parameters of the 
MSSM, by taking into account only the electroweak contributions. The growing of the 
width with tan/? (fig. ^ (a)) shows that the bottom Yukawa coupling plays a central role 
in these contributions. The evolution with the trilinear coupling A b and the higgsino mass 
parameter [i -the two parameters that appear in the trilinear coupling b^ h- displayed 
in figs. |3] (b) and (c) shows that this parameters can enhance the width some orders of 
magnitude. We have artificially let Ab grow up to large scales (that are not allowed if one 
wants that squarks do not develop vacuum expectation values) in order to emphasize the 
dependence on A b . The various spikes in these figures reflect the points where the form 
factors change sign, whereas the shaded region in fig. [3] (c) reflects the exclusion region of 
\x by present LEP bounds on the chargino mass. 

In all these figures the ratio (|3|) is smaller for the heaviest scalar Higgs {H°) because 
with the parameters @ the CP-even Higgs mixing angle a is near — vr/2, so making the 
couplings of H° with down quarks and squarks much weaker, but in fig. ^ (d) it can be 
seen that when the pseudoscalar Higgs mass grows (and this shifts a far away from — n/2) 
the two scalar Higgs change roles. 



4 



We conclude that the typical value of the ratio (0), at large tan j3 < 50 and for a SUSY 
spectrum around 200 GeV, is 



B 



SUSY-EW 



(t ch) ~ 10" 



(6) 



This is an improvement of the previous results|J, specially in the A channel, by 2 
orders of magnitude. 



4 SUSY-QCD contributions 

The gluino-mediated supersymmetric strong interactions in the MSSM can also produce 
FCNC processes. This occurs when the squark mass matrix does not diagonalize with the 
same matrix as the one for the quarks. We introduce then intergenerational mass terms 
for the squarks, but in order to prevent the number of parameters from being too large, 
we have allowed (symmetric) mixing mass terms only for the left-handed squarks. This 
simplification is often used in the MSSM, and is justified by RGE analysis^ . 



The mixing terms are introduced through the parameters 8ij defined as 



(m: 



LL)ij 



(7) 



where m; is the mass of the left-handed % squark, and m?- is the mixing mass matrix 
element between the generations i and j. These parameters are constrained by low energy 
data on FCNC [13]. The bounds have been computed using some approximations, so they 
must be taken as order of magnitude limits. We use the following bounds on the 8 



parameters [14 



\5i 2 \ < .1 ^rriu m £ /500 GeV 
|(5i 3 | < .098 ^/rriu m t -/500 GeV 
| $23 1 < 8.2mcm f /(500GeV) 2 



(8) 



We compute the contributions to diagrams of fig. [I] with gluinos (g) and up squarks (u c 
a = 1, 2, . . . , 6). We use the same input parameters as in the electroweak contributions (|5[ 
except for the n parameter, namely 



H = -200 GeV 

m-g = 150 GeV 

/ 0.03 0.03 
8 = 0.6 

V o 



(9) 



A comment is in order for the present set of inputs: we have introduced in (|5]) the 
lightest stop mass as an input, and this stop is mostly a However, in this new 
parametrization we introduce this mass as the lightest u a mass, which will be mostly a 
Ir. On the other hand, we use a larger absolute value for the /x parameter to enhance the 
ratio (|). 

Again the largest contribution comes from the right-handed form factor of (|2|), but this 
is only because we have chosen not to introduce mixing between right-handed squarks. 
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We have plotted the evolution of the ratio (0) with some parameters of the MSSM in 
fig. |4j. As can be easily guessed, the most important parameter for these contributions is 
the mixing mass parameter between the 2nd and 3rd generation of left-handed squarks, 
the less restricted one of the three (eq. (||)). In fig. f| (a) it is shown that changing <5 2 3 
by 3 orders of magnitude, the ratio (Q) can increase by 7 orders of magnitude! We can 
see in fig. ^ (b) that the fi parameter also plays an important role, like in the electroweak 
contributions (fig. ||] (c)), and for the same reasons, bringing the ratio @ up to values of 
10~ 4 . Notice that the central region of \fi\ ;$ 90 GeV is excluded by present LEP bounds 
on the chargino mass. 

The evolution with the gluino mass (fig. ^ (c)) is asymptotically quite stable, showing 
a slow decoupling. Finally in fig. £| (d) we have plotted the evolution with the pseudoscalar 
Higgs mass, it is also quite stable, until near the kinematic limit for A and H°. 

We conclude that the typical value of the SUSY-QCD contributions to @, with a 
SUSY spectrum around 200 GeV, is 

B SVSY - QCD (t^ ch) ~C(10- 5 ) , (10) 

but in favourable regions of the parameter space (i.e. large fi, or relatively light gluino) it 
can easily reach values of 10 -4 . This is 1-2 orders of magnitude larger from the previous 
estimate 0. 



5 Conclusions 

We have computed the SUSY-electroweak, Higgs, and SUSY-QCD contributions to the 
FCNC top quark decay t — > ch {h = h°,H°,A ) in the MSSM, using a mass spectrum 
motivated, but not fully restricted, by model building and Renormalization Group Equa- 
tions. 

We have found that with a SUSY mass spectrum around 200 GeV, which is well above 
present boundsp], |9|, [L(J, the different contributions to this decay are typically of the order 

B SVSY - EW (t^ch) ~ HT 7 

£SUSY-QCD (t ^ c/i) „ 1(T 5 -10- 4 . (11) 

The difference of at least two orders of magnitude between the two contributions makes 
worthless to compute the interference between the two contributions, but if the limits on 
5 2 3 (eq. (||)) improve, it should be necessary to make the full computation. 

The results ( [TTj) are an improvement of the previous results 0, specially in the A 
channel, thanks to the inclusion of the h vertex. 

It would probably be difficult that this decay can be measured either at the Tevatron, 
or at the NLC, but there exists a possibility for LHC. As an example to assess the 



discovery reach of this accelerators the FCNC top quark decays into a vector boson are [15 



LHC: B(t-*cV) >5xl0" 5 

NLC: B(t^cV) > 10 3 - 10~ 4 , (12) 

where the lack of sensitivity of NLC is due to the lower luminosity. So, if the discovery 
reach for FCNC Higgs processes are not very different from that of the gauge bosons, 
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there is a possibility to measure this decay channel at LHC even if SUSY particles are 
not seen at LEPII. 
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Figure Captions 

Fig. 1 Generic one loop Feynman Diagrams contributing to t —>■ ch. 

Fig. 2 Different form factors (|2|) for the channel t —>■ c h° as a function of tan (3, with the 
typical set of inputs of eq. @ . 

Fig. 3 Evolution of the ratio @ with (a) tan j3, (b) the trilinear coupling Ab, (c) the 
higgsino mass parameter fi, and (d) the pseudoscalar Higgs mass M A o, the rest of 
inputs are given in eq.(^|). 

Fig. 4 Evolution of the ratio (|3|) with (a) the mixing parameter between the 2nd and 3rd 
squark generations <5 2 3, (b) the higgsino mass parameter /i, (c) the gluino mass m~ g , 
and (d) the pseudoscalar Higgs mass M A o, the rest of inputs are given in eqs.(^) 
and (|). 
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